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Abstract 



Using numerical solutions of single particle dynamics, we consider a 
chaotic electron-whistler interaction mechanism for enhanced diffusion in 
phase space. This process, when applied to parameters consistent with the 
Jovian magnetosphere, is a candidate mechanism for pitch angle scattering in 
the lo torus, thus providing a source of auroral precipitating electrons. 

We initially consider the interaction between two oppositely directed 
monochromatic whistler mode waves. We generalize previous work to in- 
clude relativistic effects. The full relativistic Lorentz equations are solved 
numerically to permit application to a more extensive parameter space. We 
use this simplified case to study the underlying behaviour of the system. For 
large amplitude monochromatic waves the system is stochastic, with strong 
diffusion in phase space. 

We extend this treatment to consider two oppositely directed, broad band 
whistler wave packets. Using Voyager 1 data to give an estimate of the whistler 
wave amplitude at the To torus at Jupiter, we calculate the degree of pitch 
angle scattering as a function of electron energy and initial pitch angle. We 
show that for relatively wide wave packets, significant pitch angle diffusion 
occurs (up to ±25°), on millisecond timescales, for electrons with energies 
from a few keV up to a few hundred keV. 
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INTRODUCTION 



The electron-whistler interaction has been considered as potential mechanism for pitch 
angle scattering in planetary magnetospheres. Whistler waves are able to resonate with 
electrons over a broad energy range, from less than 100 keV to several MeV (Horne and 
Thornc, 1998). In particular the Hamiltonian has been obtained for relativistic electrons 
interacting with a whistler mode wave of single k, revealing underlying behaviour that is 
dynamically simple (Laird, 1972). 

Gyroresonance processes with near parallel propagating whister waves have been con- 
sidered as a pitch angle scattering mechanism (e.g. Kennel and Petschek (1966), Lyons 
and Williams (1984)). These processes consider electrons at resonance with many whistler 
waves of different frequencies to produce pitch angle diffusion (Gendrin, 1981). Instead we 
consider a process that is enhanced over (and in addition to) 'diffusion' associated with 
particle dynamics strictly at resonance. We consider a stochastic interaction mechanism, for 
which particles, that are not at resonance with the waves, effectively diffuse faster than the 
particles that are in resonance. For clarity, we will refer to this as 'off-resonance' diffusion. 

Stochasticity has been introduced by coupling the bounce motion of the trapped elec- 
trons with a single whistler (Faith et al, 1997), while the presence of a second, oppositely 
directed whistler wave was shown from the non-relativistic equations of motion to introduce 
stochasticity into the system, demonstrated numerically for a wave frequency of half the 
electron gyrofrcqucncy by Matsoukis et al (1998). This mechanism has been shown to exist 
in self-consistent simulations (Devine et al, 1996). Diffusion processes are also effected by 
inhomogenieties in the medium (see for example Helliwell and Katsufrakis (1974)) though in 
the case considered here variations in the background field over the short distances traveled 
by the electrons during the rapid interaction are small and the background magnetic field 
can be assumed to be uniform. 

Wave propagation calculations by Wang et al (1995) have shown that lightning induced 
whistlers are guided into the lo torus by the strong local peaks in density. The waves tend 
to oscillate around the field line corresponding to the peak in torus density (L=5.7, Wang 
et al (1995)), in a similar way to ducted whistlers with low normal angles (< 10°) at the 
equator. We consider the interaction of two oppositely directed whistlers with the same 
frequency, propagating parallel to the background field, since small differences in the wave 
frequencies and slightly oblique waves (wave normal < 10°) result in small perturbations to 
the dynamics of the system. 

In this paper we generahze the work of Matsoukis et al (1998) to consider a range of 
wave frequencies below the gyrofrcqucncy and include relativistic effects. We consider the 
efficiency of the mechanism in scattering electrons in 'pancake distributions' i.e. distributions 
with a large peak in the distribution function at 90° pitch angles (sometimes represented as 
T± 3> T||). Since we are considering single particles dynamic over a range of phase space, 
and not the evolution of the distribution function, per se, we study pancake distributions by 
considering electrons with high initial pitch angles i.e. v_i_ > v\\. These then evolve to explore 
phase space. Recent plasma density models have shown that anisotropic distributions are 
required to fit the observed whistler dispersions in the Jovian magnetosphere (Crary et al, 
1996). We investigate the dependence of the degree of stochasticity of the system (using 
Lyapunov exponents) on the wave amplitude, wave frequency and perpendicular velocity. 
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Wc simplify the equations of motion in the hmit of low wave amplitudes and extend the 
process to study two oppositely directed, broad band whistler wave packets. We obtain an 
estimate for the whistler wave amplitude using the analysis of Voyager 1 data (see Scarf et al 
(1979), Kurth et al (1985) and Hobara et al (1997))) and estimate the degree of pitch angle 
scattering for different electron energies resulting from the interaction with wave packets of 
different bandwidths. 



EQUATIONS OF MOTION 

We consider a total magnetic field of the form 

B = Bo + B+ + B- 

where Bq = -Bqx is the background magnetic field and B+ and B~ are the whistler waves 
propagating parallel and anti-parallel to the background field respectively (for coordinate 
system see Figure 1). We assume that the background field lines are uniform, since, as we 
will see, the interaction is sufficiently fast so that changes in the background field experienced 
by the electrons are small, e.g., for electrons close to Jupiter's magnetic equator at 6i?j, the 
field changes by less than 1% for an MeV electron traveling at 0.9c and interacting with the 
field for 1000 electron gyroperiods (0.1s). 
The wavefields B+ and B~ are given by 

B^=Btj[cos(A;x — ujt)y — sin(A;a; — ti;t)z] 
Bj=B(^[cos(— /ex — cut + 6'o)y — sm[—kx — ujt + 6*0)2] 

with X parallel to the background field and y and z perpendicular. The wave amphtude is 

given by B^^, while ^0 = is the initial phase difference of the waves. The wave frequency, cj, 
and wave number, k, are given by the electron whistler mode dispersion relation (neglecting 
ion effects): 

where ujpe is the plasma oscillation frequency and is the electron gyrofrequency. Electrons 
traveling at the correct parallel velocity will experience a constant field and will interact 
strongly with it. This resonance velocity, Vr = v^x is given by the resonance condition: 

a; — k • Vr = nJle/7 (2) 

where n is an integer, and 7 = (1 — v"^ /(?) ^^"^ is the relativistic factor. The corresponding 
electric field is obtained from Maxwell's relation for plane propagating waves, k'Ei^ = a;kABt^ 
and the dispersion relation (1). 

We write v — v\\'k-\-v±_ cos y + v_l sin z, where = 0(i) is the phase of the perpendic- 
ular velocity, define the phase angle ip = kx — ut + (f) as the angle between the perpendicular 
velocity and B+ and define the phase difference A-^ = ^0 as the angle between the two 
waves. 
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We substitute these into the Lorentz force law to give the full equations of motion: 



dv\\ bv± /_ (^v\\\ . , bv± /_ ujv\\\ . , , ^ 
dv 



dt 




sin(V' + AV') (4) 



dip , I b f uj\ , 

—— — KV\\ — U -\ V\\ COS ip 

dt 7 -jv± \ kj 

where b — B^^/Bq is wave amplitude scaled to the background field, and time and velocity 
have been rescaled with respect to the gyrofrequency, fie, and the phase velocity, = w/k, 
respectively. In the non-relativistic limit 7^1 these reduce to the non-relativistic equations 
of motion given in Matsoukis et al (1998), (Equations (5)-(7)). 

Reduced Equations 

We can reduce the full equations of motion (3)-(6) to demonstrate the underlying prop- 
erties of the system for low wave amplitudes. The reduced equations also give qualitative 
information on how the system will be effected when different parameters are varied. To 
simplify the full equations of motion we rescale the velocities with respect to the normalized 
wave amplitude, b: 

v\\ = vwfl + bv^i (7) 
v± = ^^±,0 + ^''^±,1 (8) 

(9) 

We substitute into Equations (3)- (6) and differentiate with respect to time. To first order 
in b, W||^o and v±fi are constant, 70 = ^1 — f | q/ — v'j_ q/ is a constant and we have 
■0 = A;x — (1/7 — a;)t, giving the following reduced equation: 

(Px ^ 2bv^ ^.^ _ _ 

dt-' 7o 

Thus we have the equation of a coupled pendulum with variables v\\ = x and x. Perturba- 
tions in and x are proportional to the wave amplitude, b, to I/70 and to v±fi. 

WAVE PACKET APPROXIMATION 

The analysis of Voyager 1 data in Scarf et al (1979) and Kurth et al (1985) shows that 
whistler mode waves at the lo torus are broad band waves. Instead of a single pair of 
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oppositely directed monochromatic whistler waves, we consider a sum of N such pairs, with 
each pair having a different frequency and (opposite) wave number. For N pairs of waves, 
the reduced equation (10) becomes: 



N 



2vi_ 



^ -hi sin [(1/70 - uji)t\ cos [hx - 60, i] 



i=l 



7o 



As N 



00 we can replace the sum by an integral over cu: 

(fx _ 2v±o, 



7o 



-b{uj) sin [(1/70 — uj)t\ cos [k(uj)x — ^o('^)] 



(11) 



(12) 



We approximate the wave amplitude using a 'top hat' distribution, i.e. we have a wave 
packet with central frequency, ujq, and constant wave amplitude across a narrow range of 
wave frequencies, Acu: 



h{u) 





ho 




UjQ-^<U<UJo + ^ 



(13) 



We approximate the dispersion relation as k{uj) = /cq + I3{uj — tfo), close to the central wave 

frequency, ujq and wave number, /cq, where (5 — 2koc'^llo'-]')°'^ ^ ^ where 1/(3 is the group 
velocity of the wave packet. Substituting in and integrating with respect to the the wave 
frequency, cu, gives the equation of motion for an electron interacting with two oppositely 
directed wave packets: 



(Px 26ow±,o . 

— = sm 

dV 7o 

, p, 2hov±fl . 

+iZe sm 

7o 



[ U!)t + kgX 

7o 



7o 



— U!)t — kox 



sin [{t - /3x)Auj/2] 
(t - (3x) 

sin [{t + 3x)X^/2\ 
{t + (3x) 



(14) 



In the limit of narrow bandwidths, Auj — > 0, the wave packet equation (14) yields the 
reduced equation (10) with wave amplitude h' — hoVleAuj. A more detailed derivation of 
the reduced equation (10), and the wave packet equation (14), can be found in Wykes et al 
(2000). 



NUMERICAL RESULTS 

Figure 2 shows numerical solutions of the full equations of motion. The phase diagrams 
are stroboscopic surfaces of section Tabor (1989) constructed from cut-planes where x — 
{n + l/2)7r/A;, to sample the full electron phase space. The initial parallel velocity, v\\, was 
varied over the range [— fr,fr], where is the resonance velocity, given by the resonance 
condition (2), to give a good coverage of phase space. 

All electrons were given the same initial perpendicular velocity, v±_ = 0.7c {v±/vr ~ 20, 
to give high initial pitch angles in order to study a pancake distribution), and the phase 
angle, ip (defined as the angle between the perpendicular velocity and the first whistler wave 
B+, see Figure 1) was initially set to or tt. 
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For low wave amplitudes, Figure 2 a), the trajectories are essentially regular and char- 
acterized by two sets of resonances. Regular trajectories are described by KAM surfaces 
(Tabor, 1989) i.e. the trajectories are near-integrable and there is an approximate constant 
of the motion associated with each orbit. The regular trajectories are represented in the 
phase diagrams as smooth lines. 

As the wave amplitude is increased in Figures 2 b) and 2 c). stochastic effects are in- 
troduced into the region between the two resonances (stochastic trajectories appear as a 
spread of dots). For higher wave amplitudes. Figure 2 d), the system becomes increasingly 
stochastic with regular trajectories confined to KAM surfaces close to the resonances. The 
dependence of the degree of stochasticity on the wave amplitude is supported by the form of 
the reduced equation (10), where the perturbations in parallel velocity scale with the wave 
amplitude. 

In Figures 2 b)-d), the stochastic regions are bounded by the first regular, untrapped, 
trajectories away from resonance. In Figures 2 c) and 2 d), there are no regular trajectories 
separating the resonances and stochastic electrons can diffuse throughout the stochastic 
region, resulting in large changes in pitch angle. It can be seen that the stochastic electrons 
are able to diffuse to a greater extent in phase space than the electrons on regular trajectories. 

In Figure 3 we show a sequence of phase diagrams for increasing perpendicular velocity, 
with constant wave amplitude, b = 0.003, and wave frequency, cu — 1/2. The reduced 
equation (10) describes pendulum like behaviour with oscillations in v\\ proportional to both 
the wave amplitude and the perpendicular velocity, v±. The resonance condition (2) shows 
that Vr = Vr{'y{v±,v\\),k,uj), therefore (for v\\ << v±) Vr, in addition to the total electron 
energy, E, can be described as a function of the perpendicular velocity when u and k are 
constant. By varying v±, only we can consider the dependence of the degree of stochasticity 
on v± and hence E. 

In Figures 3 a) and c,) the perpendicular velocity increases to relativistic velocities, v±— 
0.3-0. 75c (Energy 25-220 keV). The degree of stochasticity increases with v±. From the 
resonance condition (2) we see that increasing v± increases 7 and hence reduces Vr and the 
separation between the two resonances. 

In Figure 3 d) where v± — y/1 — — 0.86c (E=380 keV) the resonance condition (2) 
is satisfied for Vr = 0. Increasing v± further causes the resonances to pass through the 
v\\ = line and change sign. In Figure 3 e) we have v± = 0.95c (E=740 keV). The resonance 
velocity now increases with v± and the degree of stochasticity decreases, until the system is 
no longer dominated by stochastic trajectories (Figure 3 f) with v± = 0.99c, E=1.8 MeV). 

The presence of a peak in the degree of stochasticity and its dependence on v± is a 
relativistic effect. For non-relativistic velocities, 7, and hence the resonance velocity, is 
constant and the degree of stochasticity continually increases with v± (Matsoukis et al, 
1998). 

In Figures 4 and 5 we show phase diagrams similar to Figure 2 except that we now plot 
pitch angle against phase angle. The initial conditions are; b — 0.003, lo — 1/2, v± — 0.7c 
and E=175 keV. The phase diagrams are qualitatively similar to Figure 2 and share many 
of the same features. Regular trajectories are confined to close to the resonance pitch 
angle, = arctan (f_|_ o/fr), where v±^o is the initial perpendicular velocity (0.7c) and Vr 
is the resonance velocity. The stochastic region is bounded by the first regular, untrapped, 
trajectories away from resonance. Again, we see that stochastic electrons are able to diffuse 
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in pitch angle to a greater extent than resonant electrons. Later, we shall investigate the 
change in pitch angle for wave amplitudes obtained for the Jovian magnetosphere, as a 
function of electron energy, initial pitch angle and wave frequency. 

In Figure 5 diffusion in pitch angle is very fast. Significant pitch angle diffusion occurs 
on timescales of the order of tens of electron gyroperiods. On this timescale, electrons at 
Jupiter's magnetic equator (L=6), traveling parallel to the background magnetic field at 
high rclativistic velocities {v ^ c), experience changes in the background magnetic field of 
less then 1 %, therefore the approximation that the background magnetic field is uniform is 
valid. 

In Figures 6 a) and 6 b) we show sample numerical solutions of the wave packet equa- 
tion (14) for narrow wave packets {Au = Qe/500). In panel a) there is little change in pitch 
angle with time and in panel b) the trajectory is qualitatively similar to regular trajectories 
in Figure 4. In Figure 6 c), for wider wave packets, (Aa; = r2e/50), the change in pitch angle 
is greater (Aa ^ 15°, ao = 90°), with perturbations in pitch angle decreasing with time. 
The trajectory in panel d) is similar to stochastic trajectories in Figure 4. 

We have investigated solutions of the full equations of motion in order to understand 
the underlying behaviour of the system. We now consider numerical solutions of the more 
realistic wave packet equation for applications to the Jovian magnetosphere. 



ESTIMATION OF WHISTLER WAVE AMPLITUDES. 

We use the analysis of Voyager 1 data in Scarf et al (1979), Kurth et al (1985) and 
Hobara et al (1997) to estimate the whistler wave amplitude. For the lo torus at Jupiter the 
electron gyrofrequency, fig = 53.2 kHz, corresponding to a background field, Bq = 302 nT. 
At L=6 the plasma frequency, ujpe = 355 kHz (see Scarf et al (1979) and Kurth et al (1985)). 
We consider broad band whistlers with frequencies up to the electron gyrofrequency and 
bandwidths of the order of a few kHz (Ao; — ile/SO, see Hobara et al (1997)). 

The plasma wave instrument on Voyager 1 measures the electric field spectral density 
(E^/Auj) of the whistler waves over a set of frequency channels of finite width Aou. We 
estimate the corresponding magnetic field strength using: 

where = ui /k is the phase velocity given by the dispersion relation (1). Using a spectral 
density of 10"^^ V^m~^/Hz obtained by Hobara et al (1997) we obtain a wave amplitude, 
h — B^^/Bq = 10~^ for a wave frequency, uj = Qc/2. While this is insufficient for stochasticity 
for a single pair of monochromatic whistlers, for broad band whistler wave packets, significant 
diffusion does occur. 

We obtain a similar estimate using Ulysses data for the Jovian magnetopause (Spectral 
density = lO'^^'^ V^m'VHz for uj = Qe/2, see Tsurutani et al (1997)). To determine 
uniquely whether or not this process is significant, for a given bandwidth, wave amplitude 
rather than spectral density measurements are needed. 

In this context it interesting to note that for the Earth, direct amplitude measurements 
(6 = 4 X 10~^, see Nagano et al (1996)) and the extrema of spectral density measurements 
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(6 = 5 X 10~^, see Parrot and Gaye (1994)) yield whistler amplitudes sufficient for stochastic- 
ity with a pair of monochromatic whistlers, whereas average spectral density measurements 
do not (Tsurutani et al, 1989). 

Pitch Angle Diffusion in the lo torus at Jupiter 

In Figure 7 we plot the change in pitch angle for electrons during an interaction with two 
oppositely directed whistler wave packets with relatively wide bandwidths (Acu = 1 kHz, see 
Hobara et al (1997)), central wave frequency, Uq = f2e/2, and wave amplitudes consistent 
with the Jovian magnetosphcrc (6 = 10^^). Wc estimate the change in pitch angle for 
electron trajectories with initial pitch angles in the range cto = [0°, 180°] and electron energies 
from — 400 keV. Significant pitch angle diffusion occurs, from a few degrees for low initial 
pitch angles, up to 25° for higher initial pitch angles, for a range of electron energies. 

Increasing the central wave frequency has the effect of shifting the entire structure shown 
in Figure 7 down in energy. Lower wave frequencies arc required to efficiently scatter high 
energy electrons {E > MeV). Given an even distribution of waves with frequencies up to 
the electron gyrofrequency, electrons with energies from a few keV up to a few hundreds 
of keV are most readily scattered. The pitch angle diffusion mechanism described above 
is consistent with the framework described by Kennel and Petschek (1966). However, in 
their derivation of the diffusion coefficient (see their equations (3.6) — (3.10)) they make an 
estimate of the timescale for diffusion based on the dynamics of particles near resonance. 
Here, we examine the detailed dynamics in phase space and find that particles off resonance 
are chaotic and as a consequence diffuse much faster than estimated by Kennel and Petschek 
(1966). Hence we have termed this process 'off- resonance' diffusion. 

Lyapunov Exponents 

Lyapunov exponents are used to quantify the degree of stochasticity in the system. The 
Lyapunov exponents are calculated using the method described by Benettin et al (1976) and 
Parker and Chua (1994). All six Lyapunov exponents were calculated over phase space and 
evolved to their asymptotic limit. The only significant Lyapunov exponent corresponds to 
the coordinate directed along the background field. 

For positive Lyapunov exponents, two trajectories that are initially close together in 
phase space will diverge exponentially in time. For negative or zero Lyapunov exponents, 
two trajectories that are initially close together in phase space will remain close together. 
Positive (negative) Lyapunov exponents correspond to stochastic (regular) trajectories in 
phase space. For a detailed description of the theory of Lyapunov exponents see Hilborn 
(1994). 

In Figure 8 we plot the Lyapunov exponents for 10 keV electrons, averaged over a series 
of trajectories with pitch angles in the range [0°, 180°]. We consider the interaction between 

electrons and two oppositely directed whistler wave packets with relatively wide bandwidths 
{Auj = 1 kHz, see Hobara et al (1997)) and investigate the dependence of the Lyapunov 
exponent on the relativistically correct wave frequency, lu^ = jlu. The Lyapunov exponents 
are plotted as a function of l/ur^. The Lyapunov exponent is enhanced when cuj — 1/n, 
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where n is an integer. It then follows that close to these frequencies the process will be the 
most efficient in pitch angle scattering. The structure in the dependence of the Lyapunov 
exponent on the relativistically correct wave frequency is subtle and arises from higher order 
terms in the reduced equation (10), a derivation of which can be found in Wykes et al (2000). 

The separation of two electrons, with Lyapunov exponent A, in phase space, scales ap- 
proximately as ~ exp At, thus an order of magnitude estimate shows that changes in pitch 
angle scale approximately as ~ exp At, giving a characteristic time constant, r, for changes 
in pitch angle of by a factor of e^, of approximately r ~ 1/A gyroperiods. Thus the Lya- 
punov exponents in Figure 8 indicate changes in pitch angle occurring in tens of electron 
gyroperiods. 

DISCUSSION 

We have considered the interaction between two oppositely directed whistler waves and 
relativistic electrons. We initially considered monochromatic whistlers in order to under- 
stand the underlying behaviour of the system. For whistlers with wave amplitudes consistent 
with the Jovian magnetosphere the interaction with two oppositely directed monochromatic 
whistler waves results in weak diffusion in phase space. Stochastic effects can be introduced 
with high wave amplitudes: electrons not in resonance with either of the two whistlers can 
diffuse extensively in phase space. We refer to this process as 'off-resonance' diffusion as 
only the electrons not in resonance are scattered in this way. Resonance diffusion remains 
unchanged. 

We have extended the treatment of monochromatic whistlers to consider two oppositely 
directed whistler wave packets. For relatively broad band waves with amplitudes consistent 
with quiet times in the Jovian magnetosphere, we have shown that the 'off-resonance' dif- 
fusion results in significant pitch angle diffusion (up to 25°) on timescales of the order of a 
few tens of electron gyroperiods (milliseconds). 'Off-resonance' diffusion is most effective in 
scattering electrons with energies from a few keV up to a few hundred keV (MeV electrons 
require extremely low wave frequencies for effective scattering). 

We have shown that the Lyapunov exponent is enhanced when the relativistically correct 
wave frequency a;-^ = 70; = 1/n where n = 1,2,3, ... This phenomenon arises purely from 
the interaction of the two whistler waves. We expect enhanced phase space diffusion close 
to these frequencies. 

Acknowledgements W J Wykes and S C Chapman are funded by PPARC. 
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FIGURES 



FIG. 1. Illustration of the coordinate system used in the model. In a) the magnetic wavefields 
and Bj lie in the y, z plane, perpendicular to the background field, Bq = -Bqx. The phase 
angle ip is defined as the angle between B+ and the electron perpendicular velocity, v_|_ and A^/; is 
the angle between B+ and Bj. In b) the electron pitch angle a is defined as the angle between 
the velocity vector v and the background field Bq. 

FIG. 2. Stroboscopic surface of section plots for w/Og = 1/2, initial v± = 0.7c. The parallel 
velocity, v\\, has been scaled to the phase velocity, Vphase- For this wave frequency, the dispersion 

relation (1), gives Vphase = 0.07c. The phase angle ij) is defined as the angle between the perpendic- 
ular velocity and the whistler wave propagating in the +ve x direction. For low wave amplitudes, 
panel a), all trajectories are regular and the equations of motion reduce to a pendulum equation, 
with resonances given by the resonance condition (2). For slightly higher wave amplitudes, pan- 
els b) and c), stochastic effects appear as regular trajectories are broken down. For high wave 
amplitudes, panel d), phase space is dominated by stochastic trajectories with regular trajectories 
confined to KAM surfaces close to the resonances. The stochastic region is bounded above and 
below by the first regular, untrapped, trajectories away from resonance, therefore there is a limit 
on the diffusion of electrons in phase space. 

FIG. 3. Stroboscopic surface of section plots for u/Vle = 1/2 and h = B^/Bq = 0.003. Electron 
energy increases as a function v±. In panels a) - c) the degree of stochasticity increases with 
increasing perpendicular velocity (and electron energy, E) with constant co,i^e, and k. As vj_ and 
7 increases the resonance velocity, Vr, decreases (see the resonance condition (2)). In panel d) the 
resonance velocity is zero. Here v± = 0.86c and E = 300 keV. Increasing v± (and E) further in 
panel e) causes the resonance velocities to cross the = line and the degree of stochasticity now 
decreases, until in panel f) with v± = 0.99c {E = l.SMeV) the trajectories become regular once 
more. 

FIG. 4. Pitch angle plotted against phase angle for b = 0.003, uj/0,e = 1/2, v_i_ = 0.7 and 
E = 175 keV. Phase space is divided into stochastic and regular regions in a similar way to Figure 
2. Electrons with regular trajectories close to the velocity resonances in Figure 2 are confined 
to close to the resonance pitch angle = arctan(uxo/f^r)- Pitch angle diffusion is limited by 
the extent of the stochastic region, which is bounded by the first untrapped orbit away from the 
resonances. The electrons that diffuse to the highest U|| in Figure 2 have the greatest change in 
pitch angles. 

FIG. 5. Series of phase plots showing diffusion in pitch and phase angle for a small set of initial 
conditions. Diffusion is rapid with electrons reaching the minimum pitch angle in tens of electron 
gyroperiods. After 100 electron gyroperiods, electrons have diffused throughout the stochastic 
region. Parameters are as in Figure 4. 



FIG. 6. Changes in pitch angle for sample solutions of the wave packet equation for two widths 
of wave packet. In panels a) and b) Au; = Oe/5000 = lOHz resulting in very narrow wave 
packets. Trajectories are regular and resulting changes in pitch angle are low. In panels c) and d) 
Au = rie/50 = IkHz, with initial pitch angle a{t = 0) = 90°. The trajectories are now stochastic, 
with greater change in pitch angle. In panel c) it can be seen that perturbations in the parallel 
velocity decrease with time. 

FIG. 7. Pitch angle diffusion for relatively wide, oppositely directed whistler wave packets 
{Alo = r^e/SO) with central frequency to = and low wave amplitude (6 = 10~^). We plot 

the change in pitch angle against initial pitch angle and electron energy. Significant pitch angle 
diffusion occurs (from a few degrees up to 25°) for a range of electron energies and initial pitch 
angles. 

FIG. 8. We investigate the dependence of the Lyapunov exponent on the relativistically correct 
wave frequency = juj, for relatively wide band whistler wave packets. We plot the Lyapunov 
exponent as a function of The Lyapunov exponent is enhanced when lo^ = 1/n, where n is 

an integer. 



SOS : x=riiv5k , ^J^f=0 .00 1 , (i}= 1/2 , v^=0.7c SOS : x^mJ2 k, BJB^=0.00^5, ai=V2, .7c 




-+-2 024 -+-2 024 

P hasB Ang b, y Phase Ang b, y 





4 



Ph^eAngb, >|J P hasB Ang b, >|* Phase Angle, y 



SOS: x^mf^2k, BJBjj=0.003, (i^=^f2, v^=0.7c 




Time=On 



Time=20 D. 



Tir™=4CI n 



SO 



30 



80 




100 



so 



2 

Ph 



-2 
Angb, >{f 



100 



TimB=60 Q. 



2 0-2 
Phase Angb.ll* 

TimB=100 Q. 



SO 



100 



SO 



2 0-2 
PhasB Angle, l|J 

Time=1000 SI 




81.1 



81 A 



go 



.35 



□J 



I SO 



75 



Mjzzil /5000=l0hlz 




1000 2000 3000 4O00 5000 
Time 

il(ij=il /50=1 KHz 



I I 



_i I I i_ 



1000 2000 3000 4O00 5000 
Time 



37 .4 



il(ij=il /5000=lOHz 




-2 2 

Phse, \{f 

il(ij=Q /50=1 KHz 




